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Abstract. - We propose that the spin glass phase of cuprates is due to the proliferation 
of topological defects of a spiral distortion of the antiferromagnet order. Our theory explains 
straightforwardly the simultaneous existence of short range incommensurate magnetic correla- 
tions and complete a-b symmetry breaking in this phase. We show via a renormalization group 
calculation that the coUinear 0(3)/0(2) symmetry is unstable towards the formation of local 
non-coUinear correlations. A critical disorder strength is identified beyond which topological 
defects proliferate already at zero temperature. 



Introduction. - Until recently, the spin glass (SG) phase of La2-xSrxCu04 (LSCO), 
with 0.02 < a; < 0.05, was believed to have short ranged commensurate correlations, with 
incommensurate (IC) correlations appearing only at the onset of superconductivity x ^ 0.05 

New data on LSCO have however revealed IC spin correlations also in the SG regime 
The short range nature of the correlations may point to the existence of a strongly disordered 
stripe glass, in accord with our previous analysis which predicted an instability of the stripe 
array towards a disorder dominated phase. Within the picture of disordered stripes, one is 
however immediately confronted with a problem of two fundamentally different scales: the 
inter-stripe distance, £s, which scales as ^ Q and the average separation, id, between 
disorder centers (Sr ions) which scales as id ~ l/y/x. A conflict arises because for pinned 
stripes, both scales are associated with the charge sector. As id ^ is in the SG regime, the 
stripes must break up into short segments or drastically reduce the on-stripe charge density 
to take advantage of the disorder potential |6| . 

Such a conflict does not arise in models which assume that holes localize individually near 
their dopants, as indicated by transport measurements [Q, and act as frustration centers to 
the AF order [p|- [ll[|. A polarization of the frustration centers can then lead to the appearance 
of IC correlations [[12[. In such a scenario, the scale is exists only in the spin sector whereas 
id is characteristic of the charge distribution. In the following we revisit this model, assuming 
that holes produce dipolar frustration, and show that it can well explain the experiments. 
Our analysis sheds light on the origin of the cluster spin glass phase which in our analysis 
arises from topological defects of a spiral spin texture. 
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Because dipolar frustration produces a long range distortion of the AF order, this problem 
can be treated within a continuum field theory approach. A non-linear tr-model (NLctM) for 



where ps is the spin stiffness, T the temperature and n a three component unit vector rep- 
resenting the local staggered moment (we set = !)• The quenched field represents the 
dipoles, with zero mean and a Gaussian distribution. An analysis of (^) leads to a finite corre- 
lation length at T = for any finite disorder strength and to fair agreement with the tempera- 
ture and doping dependence of the correlation length of LSCO As it stands, the model 
however predicts commensurate AF correlations, in contradiction with experiments |^-^. 

It is well known that the dipole model can lead to IC correlations The Hamiltonian 
(|l|) favors the formation of a spiral phase, with a non-zero average twist of the AF order 
and a simultaneous alignment of the dipoles, (f^j) 7^ 0, as long as the lattice and spin degrees 
of freedom of dipoles are annealed. While the spiral phase has been found to be unstable 
in the case of an annealed hole distribution towards a local enhancement of the spiral pitch, 
this instability arises from the fermionic susceptibility an hence is absent in the model 
considered here, where all fermionic degrees of freedom are quenched. Further, as the lattice 
part of the vector should reflect the symmetries of the underlying lattice, a discrete set 
of favored lattice vectors for the formation of the spiral exists. Thus, the a-b (or square 
lattice) symmetry breaking associated with the formation of spiral correlations can have truly 
long range order. The continuous symmetry of spin space on the other hand inhibits long 
range magnetic order in the 2D system for either finite temperature or finite disorder. The 
experimental observation of a macroscopic a-b asymmetry [|| but very short spin correlation 
lengths thus clearly motivate the study of the dipole model. 

We here investigate the dipole model assuming a random spatial distribution of the dipoles 
but allowing for non-zero ordered moments. The non-coUinearity of the ground state which 
arises from ordered moments requires a serious modification of previous coUinear approaches 
for two reasons: First, the symmetry breaking scheme of non-collinear magnets leads to the 
appearance of three rather than two Goldstone modes which couple to disorder. Second, a 
spiral phase allows for pointlike topological defects (TD's) with Z2 charges which orig- 
inate from a chiral degeneracy of the spiral ground state and which are absent in coUinear 
models. Below, we investigate the influence of disorder on both the Goldstone modes, using a 
renormalization group (RG) approach, and on the creation of TD's, using a free energy argu- 
ment. A strongly disordered regime, as associated with a SG, is only found once topological 
defects of the spin texture are accounted for. The situation encountered here thus resembles 
that of disordered planar XY models, where the coupling of the disorder to spin waves was 
also found to lead to a simple renormalization of the spin stiffness whereas a phase transition 
to a disordered phase is driven by a proliferation of TD's |l5| . There are however profound 
differences between these simpler XY models and the present model. Specifically, the renor- 
malization of the spin stiffness in the later case is primarily driven by the curvature of the 
order parameter space while in the former case the renormalization arises from the presence 
of TD's. Further, the influence of TD's on the renormalization of non-collinear Heisenberg 
models is at present not understood and a unified renormalization group approach which ac- 
counts for both TD's and curvature terms is not available [Q. We therefore derive the RG 
equations without incorporating TD's while we address the importance of TD's in a separate 
analysis below, where we determine the critical disorder threshold of the unbinding transition. 
Below this threshold, but not above, the RG decribes the system correctly. 



this problem reads |P,[Tl|] 




(1) 
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Derivation of the model. - As in presence of spiral correlations the 0(3) symmetry of 
the spins is broken completely, the 0(3)/0(2) AF model (j^) is not adequate. We need a 
formulation of the problem which incorporates the order parameter for a spiral, which is an 
element of 0(3) x 0(2)/0(2) Thus, there are three instead of the two Goldstone modes of 
collinear magnets. A possible representation of the local spiral order is in terms of orthonormal 
nk, k = 1,2, 3, with = Skq- A derivation of a continuum field theory for a spiral state from 
a lattice Heisenberg model can be found in Using Sy — rii cos(k5 • r^j) — n2 sin(k5 • r^j) 
and n3 = rii X n2, where is the spin at the lattice site (i, j), kg = (tt, tt) + qs and is 
the IC ordering wave vector of the spiral, the effective classical Hamiltonian can be written 
in the form 

H = ^Jd^x Pk{d^,nkf + s^J d'x m • a^na (2) 

where we ignore small anisotropics of order q| in the stiffnesses pk- Roughly, pi^2 — J/i'^T) 
and p3 ~ 0. The vector s is to lowest order given by s = Jqs/T. The second term makes this 
Hamiltonian unstable, which simply shows that the pure Heisenberg model does not support 
a spiral ground state. We show below that this term is cancelled by a term originating from 
the coupling of the spins to the dipolcs. Thus, the ordered dipoles stabilize the spiral phase. 

We use a phcnomenological form for the coupling of the spiral order parameter to the 
dipoles. The dipoles locally cant the spin order and couple to the first spatial derivatives of 
the rifc fields. To generate the frustration produced by the dipoles we thus introduce a minimal 
coupling in the first term of (^), i.e. we replace (S^n^)^ with [{d^ — iB^ • L)nfe]^ where 

is a random dipole field. The components of L are 3x3 matrix representations of angular 
momenta which generate rotations about the three spin axes, with — iB^ • L nj, = B^ x n^. 
Let us write B^ = (B^)^ + so that (Q/j)^, = 0, where {■ ■ ■) jj is the disorder average. We 
obtain the following Hamiltonian for the spiral in presence of disorder. 



-f^ = ^ J d^xpkid^Uk)'^ + J Pkd^nu ■ x rifc. 



(3) 



where the ordered part of the dipole field cancels the second term in (Q). Thus, pkd^rik ■ 
(Bfj) X rifc + Spiii ■8^^112 = 0. As qs oc s, this equation relates the incommensurability linearly 
to the density of ordered dipoles. The remaining part of the dipole field, Q^, is a quenched 
variable with zero mean and we assume Gaussian short ranged statistics, {Qfi{^)Qt{y)) = 
A^(x — y) Sab S^u. In absence of disorder, the Hamiltonian is 0(3) x 0(2)/0(2) symmetric. 
The model reduces to in the case pi,2 = with p^ = Ps/T, — n and ~ x n. 
In Fourier space, the disorder coupling takes the form of a correlated random field coupling 
/ nfc(-q) • hfc(q), with 

hfc(q) = ipkq^, j d^x (Q^ x nfe)e*'i '' (4) 



and correlations ^/i^(q)/i^, (q')y oc ^(q— q')|q| • Using an Imry-Ma type argument it can be 

shown |pj| that this coupling is marginal in two dimensions and that therefore a RG analysis 
is required to investigate the scaling of the disorder term. 

Renormalization. - The RG is obtained in the SU{2) representation p9[ | where 7iJ5 = 
trja°-ga''g^^] /2 (cr" are Pauh matrices and g e SU{2)). Introducing = -itr[a°-g^^dfj^g] /2 
pO| , eq. (|^) can be written in the form, 

H=jjd^x [aI + aAf] + 2jd^x Pk e,jk fabc n"^ nl (5) 
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where t ^ = 2{pi + pa) and a = (pi —pz)/{pi +^3)- At the point a — Q the symmetry is 
enhanced to 0(3) x 0(3)/0(3) ~ 0(4)/0(3) while at a = -1 the model is collinear. The one- 
loop RG equations are obtained by splitting the field g into slow and fast modes, g — g' exp(i 
(jf'a"') and tracing out the fast modes (/)° with fluctuations in the range [A, 1] (where the 
original UV cutoff is 1). For A <C i the RG equations are given by 

AT:dtt = 2{l - a)t^ + [2 - a + a^)Xt (6) 
Airdga = -4a(l + a)< - a(l + a)(3 - a)A 

where £ = In (A). For A = 0, these equations describe the RG of a clean spiral while 
for the collinear point a = — 1, the equations reproduce the RG of the stiffness for disordered 
collinear models From (||) it is seen that there are two fixed points for a, see also fig. 0. 
The collinear point a = — 1 is unstable whereas a = is stable, irrespective of the disorder. 
The coupling to weak disorder thus only renormalizes the spin stiffness but does not not lead 
to any new fixed points. 

The validity of is however limited to high temperatures, for low temperatures the 
RG of A must be accounted for. We have calculated the RG of A using the method described 
in However, for a ^ new disorder coupling terms appear. Specifically, the RG generates 
transverse correlated fields which couple to the ni_2 vectors even at a = — 1 (i. e. pi^2 = 0), 
whereas in the original Hamiltonian with a — —1 they only couple to na, see eq. (^). Such 
transverse correlated fields destroy the collinear fixed point. Thus, even if the original AF 
order is collinear (i.e. in absence of dipole ordering), the disorder drives the system to a 
non-coUinear state. We recover the RG equations for the collinear model obtained in only 
if we ignore the correlated transverse field coupling. A purely collinear analysis is thus not 
valid in presence of dipoles and cannot describe the low temperature regime correctly. This 
is in fact expected, because a random canting of the spins destroys the remaining 0{2) spin 
symmetry of the collinear model. We conjecture that the system will always flow to the point 
a — 0, although this remains to be proven. 

At the highest symmetry point a = however, no new coupling terms are generated. It is 
then straightforward to arrive at the following RG equations, valid for a = but any initial 
ratio of A/t, 

2TTdit = t^ + \t ; indiX = A^ (7) 

which can be simplified through z = t + X/2 to get 2iTdiz — z"^. So for a — Q the presence 
of disorder leads to an additive renormalization of the stiffness, t ^ t + X/2. In presence of 
any amount of disorder, the IC correlation length ^ at T = is finite, as can be inferred from 
an integration of the RG equation, yielding ^ cx exp(CAo ^) at T = with some constant C. 
While the disorder scales to strong coupling, the relative disorder strength with respect to 
the stiffness, X/t, always scales to zero so that at long wavelengths the disorder becomes less 
relevant. This is surprisingly different to the situation with a = — 1 fixed |jll[], where the ratio 
X/t was found to diverge below a certain initial value of Ao/to which was interpreted as the 
scaling towards a new disorder dominated regime. Thus, if one correctly takes into account 
non-coUinearity, this cross over to a strongly disordered phase disappears. 

To summarize, within the perturbative RG analysis disorder leads only to a simple renor- 
malization of the spin stiffness which in turn leads to a finite correlation length already at 
T = 0. As mentioned already in the introduction, disorder may have a more dramatic effect 
via the creation of TD's in the spin texture. Below, we analyse the stability of the spin system 
against the creation of such defects. 
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Fig. 1 - RG flow for A < 



Fig. 2 



Fig. 2 - Single topological defect of a spiral with a > 0. The arrows represent the local staggered 
moments and the spin vector of ordered dipoles point out of the plane. 



Topological defects. - TD's in spirals are related to the chiral degeneracy of the spiral, 
i.e. a spiral can turn clock- or anti-clock wise These defects are pointlike and at a TD 

the spiral changes its chirality. Saddle point solutions of TD's for spirals are of the form |^ 

5(x) = exp[imV*(x)/2], (8) 

where m is a space independent unit vector and ^'(x) satisfies V^5'(x) = with singular 
solutions ^'(a;,?;) — arctan(j//a;). In fig. |^ the spin distribution around an isolated defect 
in a homogeneous spiral is shown for a > 0. The energy of a TD diverges logarithmically 
with the linear system size R, (3E = [1 + (rn^)^a]7rlni?/2f, so for a < the lowest energy 
defects have = ±1 while for a > 0, to^ = is preferred. Because of the logarithmic 
divergence of the energy, isolated TD's are absent in clean systems at sufficiently low enough 
temperatures and in this regime the predictions of the NLcrM analysis should hold. Indeed, 
NLctM predictions for related frustrated Heisenberg models were found to be in excellent 
agreement with numerical simulations at low temperatures [^3| . Only at higher temperatures 
an abrupt decrease of the correlation length was observed which has been attributed to the 
appearance of unpaired TD's [^ , |23|] . It has been argued, that the high-temperature regime 
is similar to the equivalent phase of XY models, as both are governed by the plasma phase 
of the Coulomb gas model . The spiral defects differ however from their XY counterparts 
in that the logarithmic divergence of the energy appears only for length scales smaller than 
the correlation length obtained from the NLctM analysis. In an expansion around the saddle 
points the defects do not decouple from the spin waves even at the lowest (Gaussian) level. 
Also, TD's of spirals carry a Z2 charge, whereas XY vortices have Z charges. Despite these 
differences, the critical temperature, at which free TD's first appear can in both cases be well 
estimated using a free energy argument p^ . 

A comparison to XY models is further usefuU as the disorder coupling we employed is very 
similar to the one appearing in XY models with randomly fluctuating gauge fields In 
these models, if one ignores vortices, the influence of the disorder was shown to amount to a 



simple renormalization of the spin stiffness |15|, and no disordering transition as a function 
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of the disorder strength is found. Once TD's are inchided in the analysis, randomness leads 
to a disordered phase even at T = through the creation of unpaired defects if the disorder 
fluctuations are stronger than some critical value (l^,^ . Similar to the estimate of the critical 
temperature at which thermal fluctuations lead to the appearance of TD's , the critical disorder 
strength is well estimated from the free energy of an isolated defect in presence of disorder 
26|. Here, we use this approach to calculate the critical disorder strength for the spiral 



state. The free energy of a defect state is given by PF = (1 + (m^)^a)7r In i?/(2t) — (InZ^;^, 
where the second term contains the contribution of disorder, 

Zd = J (fyexp J Pk eijk^abc A'-^n'^nlQ''^ 

with A^, rifc derived from (H). Using the replica trick one finds at low temperatures I3F = 
2pi7r(l - y/SX/n) \nR for a < and f3F = {pi + P3)n{l - ^SA/tt) Ini? for a > 0. Thus, at 
low temperatures and A < Ac = 7r/8 free defects are absent and the system is well described 
within the NLcrM analysis above. For strong disorder, A > Ac, the creation of TD's is favorable 
even at T = 0, leading to a proliferation of free TD's and much shorter correlation lengths 
than predicted from the NLcrM, without affecting the incommensurability. Anticipating the 
analogy to XY models, for A > Ac we expect an IC correlation length roughly given by 
^ oc exp(&/v^A — Ac) (with some constant b) [^. However, instead of the divergence of ^ at 
A = Ac a crossover to the NLcrM result is expected. 

Comparison with experiment. - Experiments found the incommensurability to scale 
roughly linearly with doping which points to a doping independent ordered fraction of dipoles. 
To judge, whether or not TD's play a role in the LSCO SG phase, we must estimate A. We 



can use as a lower bound for A the result obtained from the coUinear analysis |11| where a 
disorder parameter equivalent to ours, but defined on the much smaller scale of the AF unit 
cell, was used. From a fit of the correlation lengths behavior at x < 0.02 and large tempera- 
tures, it is obtained A ~ 20a;. In this regime of x, the low temperature phase has long range 
AF order and a coUinear analysis is well justified. This relation should remain valid also in 
the SG regime, a view consistent with experiments where the width of the distribution of 
internal magnetic fields (i.e. local staggered moments) was found to increase simply linearly 
with doping, with no detectable change on crossing the AF/SG phase boundary If we use 
our above estimate for the critical disorder strength we find remarkably a critical Xc ^ 0.02. 
Considering that A ~ 20a; is a conservative lower bound of A at the long length scales relevant 
to spirals, we conclude that in the entire SG phase, free TD's will be present already at T = 0. 
Thus, we expect at T = IC correlation lengths which roughly follow ^ exp(6/V20a; — Ac). 

Our suggestion, that the IC correlations are related to ordered dipolar frustration centers 
can be directly tested experimentally in co-doped samples La2_xSrxCui_zZn204. It was 
argued that co-doping has the effect of renormalizing the density of dipoles to lower 
values X x{l — jz) (experiments indicate 7 ~ 2 |Q) as Zn atoms placed close to a locafized 
hole strongly distort the hole wavefunction which leads to a reduction or complete destruction 
of the AF frustration induced by the hole. Thus, co-doping with Zn has two effects: First, it 
lowers the amount of frustration in the sample and thus enhances the correlation length |^ , |29[| . 
Secondly, Zn doping also lowers the total amount of ordered dipoles which leads to a decrease 
of the incommensurability by a factor 1 — 72 which should be observable experimentally. 

Conclusions. - We propose a novel description of the SG phase as a strongly disordered 
spiral state, which can account for the IC correlations and the a-b asymmetry without invoking 
any kind of charge order. We show that a coUinear analysis is inadequate to describe the 
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strongly disordered regime. No sharp transition towards a disorder dominated phase is found 
within the NLcrM model analysis and TD's of the spiral must be accounted for to explain the 
extremely short correlations observed in experiments [ ^0| . The picture of the cluster SG phase 
as proposed in JlO| has some similarities to our work in that it invokes the presence of defects. 
However, in this picture, the defects are those of an 0(3)/0(2) spin system (skyrmions). Also, 
while the skyrmion model predicts IC correlations only for x > 0.05 ]lOt , our analysis includes 
IC correlations from the outset. 
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